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Abstract. We model the nucleosynthesis during a radiative interpulse phase 
of a rotating 3M Q Asymptotic Giant Branch (AGB) star. We find an enhanced 
production of the neutron source species 13 C compared to non-rotating models 
due to shear mixing of protons and 12 C at the core-envelope interface. We 
estimate that the resulting total production of heavy elements by slow neutron 
capture (s-process) is too low to account for most observations. This due to 
the fact that rotationally induced mixing during the interpulse phase causes a 
pollution of the 13 C pocket layer with the neutron poison 14 N. As a result we 
find a maximum neutron exposure of r ma x = 0.04mbarn _1 in the s-process layer 
of our solar metallicity model with rotation. This is about a factor of 5 ... 10 less 
than required to reproduce the observed stellar s-process abundance patterns. 
We compare our results with models that include hydrodynamic overshooting 
mixing, and with simple parametric models including the combined effects of 
overshooting and mixing in the interpulse. Within the parametric model a range 
of mixing efficiencies during the interpulse phase correlates with a spread in the 
s-process-cfficicncy. Such a spread is observed in AGB and post- AGB stars as 
well as in pre-solar SiC grains. 



1. Introduction 

The s-process is characterized by neutron captures that are immediately followed 
by /3-decays. Unstable isotopes involved in the s-process path have typical life 
times of the order of hours; at neutron densities of N n < 10 10 cm -3 they decay 
rather than capture another neutron (see Busso etal., 1999, for a general review 
of the s-process). Current models (Gallino etal., 1998; Goriely & Mowlavi, 2000; 
Lugaro etal., 2002) favour low mass Asymptotic Giant Branch (AGB) stars 
as the dominant production site for the s-process nuclei with mass numbers 
A > 90. In evolved AGB stars He-shell burning is periodically unstable and 
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Figure 1. s-process indices [s/Fe] and [hs/ls] as a function of neutron ex- 
posure. As explained in the text this figure translates observational properties 
of s-process enriched stars into properties of the partial mixing zone. The 
dashed/dotted horizontal lines are discussed in the text. 

generates thermal pulse (TP) cycles. The main neutron producing reaction is 
13 C(a, n) 16 at temperatures as low as T = 9 • 10 7 K. The neutrons are released 
under radiative conditions during the interpulse phase. 

The necessary amount of 13 C is produced if protons and 12 C can partially 
mix at the core-envelope interface. This may happen after TP during the dredge- 
up phase during which a fraction of intershell material is carried into the en- 
velope. The required extra- mixing is still uncertain. With the H/ 12 C-ratio de- 
creasing continuously from a few hundred in the envelope to zero in the intershell 
both a 13 C as well as a 1 N pocket forms. Herwig (2000) has described models 
with hydro dynamical, diffusive overshooting that leads to a partial mixing zone. 
However, stars rotate and in this paper we discuss the impact of rotationally in- 
duced mixing in AGB star models on the s-process. A more detailed discussion 
is available in Herwig etal. (2003). 

2. Properties of the Partial Mixing Zone 

The s-process abundances observed in AGB stars can be characterized by aver- 
aged abundance indicators. The index [hs/ls] describes the distribution of the 
s-process elements. Stars of solar metallicity show —0.5 < [hs/ls] < 0.0, where 
[ls/Fe]=^([Y/Fe] + [Zr/Fe]) and hs is the average of Ba, La, Nd and Sm. The to- 
tal overproduction of s-process elements correlates with the index [s/Fe], which 
is the average of Y and Nd, and observationally < [s/Fe] < 1 (Busso etal., 



The observed overproduction factors in the envelope are related to the over- 
production factors in the s-process zone by a number of dilution factors, which 
result from the subsequent mixing events: the He-flash convective mixing and 
the third dredge-up. The envelope abundance of any species in the envelope 
after dredge-up events in m identical TP cycles is related to the abundance in 
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the s-process zone (partial mixing zone, PM) by 

M pm Mdup /-.x 

Fenv = gmFpM (1) 

Mi S M env 

where Mpm, Mdtjp> Mis and M env are the masses of the partial mixing zone, 
the dredged up layer, the intershell zone covered by the He-flash convection and 
the envelope mass respectively. These quantities change with the pulse number, 
typical numbers in low mass stars are Mdtjp ~ 3 • 10~ 3 M Q , Mis ~ 10 -2 M Q 
and M env ~ 0.5 Mq. The factor q describes the effect of the overlapping of the 
He-flash convection zone in subsequent thermal pulses and can be estimated to 
be q = 2.3 (see Herwig etal., 2003 for details). We may assume that not more 
than 20 thermal pulses with dredge-up contribute to the envelope enrichment 
with s-process material (m = 20). 

By evaluating Eq. Q for the numbers specified above we can derive a loga- 
rithmic expression that relates the average s-process overabundance in the PM 
zone with the mass of that zone: 



logMpM = -[s/Fe] PM + c (2) 

where Y env in Eq. Q is given by using [s/Fe] en v = 1, and c = —0.44 for m = 20. 

In Fig. n w e show the variation of [s/Fe] PM and [1is/1s]pm with the neutron 
exposure in a s-process one-zone model. As more neutrons are released the 
average s-process overabundance [s/Fe] PM increases. The partial mixing zone 
should not extend into the He-shell. This puts a lower limit on [s/Fe]pM and 
thus on the neutron exposure in the partial mixing zone. In Fig. ^ the dashed 
horizontal line corresponds to Mpm < 10 -2 M Q and c = —0.44, and according to 
Eq. (J2J) [s/Fe]pM above that line are allowed. From [hs/ls] b s < 0.0 (dotted line) 
it follows that the neutron exposure should not exceed ~ 0.45mbarn _1 in the 
partial mixing zone. Otherwise the envelope s-process abundance distribution 
will eventually be top-heavy as well, and not in agreement with observations 
(Lugaro etal., 2003). 

Thus, we estimate that the neutron exposure at the end of the interpulse 
phase should be in the range 0.25mbarn _1 < r < 0.45mbarn _1 . Using Eq. (|2"|) 
and [s/Fe] PM = 4.0 at the upper r limit, the mass of the partial mixing zone 
should obey log AfpM > —4.44 (assuming m = 20). These estimates are roughly 
in agreement with the much more detailed analysis of the partial mixing zone 
by Gallino etal. (1998). 



3. s-Process nucleosynthesis in a rotating AGB model star 

Rotation induces mixing in stars through a variety of processes (Maeder & 
Meynet, 2000). In rotating AGB stars a large angular velocity gradient devel- 
ops at the interface of the compact and fast rotating core and the slow rotating 
envelope. For our analysis we take the corresponding mixing efficiencies and 
thermodynamic input from a 3Mq stellar model sequence including rotation 
(Langer etal., 1999). These models predict that shear mixing is continuously 
present at the core-envelope interface throughout the interpulse phase. For the 
s-process the effective Lagrangian mixing efficiency is important. Note, that 
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Figure 2. Abundance profiles in the partial mixing zone immediately after 
the end of the dredge-up phase (upper panel) and when the 13 C(a,n) 16 
reaction is activated (lower panel). 

the partial mixing layer is contracting significantly during the interpulse phase. 
Around the time of 13 C(a,n) 16 ignition (interpulse phase <\> ~ 0.3), the dif- 
fusion coefficient is D rot ~ 10~ 13 MQ/yr, which gives for a mixing timescale of 
~ lOOOyr a mixing range of ~ 1O- 5 M . 

We use our nuclear network and mixing post-processing code to evaluate 
the effect of mixing on the s-process. This code solves simultaneously for mix- 
ing processes and all considered nuclear processes (charged particle reactions, 
neutron captures and /3-decays). In the network all relevant charged particle 
reactions from hydrogen to 31 P have been considered with rates mainly taken 
from Angulo etal. (1999). Neutron captures for 15 light isotopes and 9 iron 
group species have been considered and most cross sections are from Bao etal. 
(2000). The neutron capture on species beyond 60 Ni are considered by a two 
step heavy neutron sink. We introduce two artificial particles ( 63 G and 1 L) with 
two extra reactions ( 62 Ni(n, 7) 63 G and 63 G(n, 1 L) 63 G) which together give an 
estimate on s-process efficiency. For the Maxwellian averaged neutron capture 
cross section of 63 G we take 0"63Q(8keV) = 120mbarn, which has been determined 
from extensive one-zone model tests with a full s-process network. 

In Fig. |21 we show the abundance profiles around the core boundary at two 
times of the post-processing calculation of one interpulse phase of the model 
sequence with rotation. The top panel shows the initial condition where rota- 
tionally induced mixing has caused some diffusion of protons into the 12 C rich 
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Figure 3. Final neutron exposure at the end of the interpulse for different 
assumptions on mixing. All cases start with a H/ 12 C partial mixing zone from 
hydrodynamic overshooting. During the interpulse no mixing (dotted line) 
and constant Lagrangian mixing coefficients have been applied (D m in units 
of M Q /yr). 



core at the end of the (in this case shallow) third dredge-up. Goriely & Mowlavi 
(2000) have shown that in such a situation of a continuously decreasing proton 
abundance profile into the 12 C rich core the mass range with a proton abun- 
dance of —2 < logXp < —3 will eventually host s-process enriched material. 
According to this criterion, and given the required properties of the partial mix- 
ing zone as estimated in the previous section the partial mixing zone obtained 
in this rotating model is too thin. 

The lower panel of Fig.|2] shows the situation when neutrons are begin- 
ning to be released by the 13 C(a, n) 16 reaction. In models without rotation, 
where the partial mixing zone results for example from hydrodynamic over- 
shooting, the 14 N and the 13 C pocket keep mainly separated. With rotation 
instead the 14 N abundance exceeds that of 13 C in the entire partial mixing zone. 
14 N is an efficient neutron absorber, most notably through the (n,p) reaction. 
The recycling effect of the protons released in this reaction and the subsequent 
12 C(p,7) 13 C reaction does not balance the poison effect of 14 N. As a result 
the neutron densities are low and the final maximum neutron exposure is much 
lower (r ~ 0.04mbarn~ 1 )than our estimates from the previous section require. 



4. Parametric Models and Discussion 

It seems that current models of rotation in AGB stars predict interpulse mixing 
coefficients that could inhibit the s-process during the interpulse phase due to 
dilution of the 13 C pocket and admixture of the neutron poison 14 N. In order to 
further study the effect of mixing during the interpulse we have constructed some 
parametric s-process models assuming a constant diffusive mixing coefficient D m 
in Lagrangian coordinates during the interpulse phase. We start in all cases 
with abundance profiles at the core-envelope interface shaped by hydrodynamic 
overshooting as described in Herwig (2000). The final neutron exposure at the 
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end of the interpulse phase is shown in Fig. for several cases. The case without 
mixing shows a larger neutron exposure than allowed according to the estimates 
in For the cases with constant D m we find an anti-correlation of mixing 
efficiency and final neutron exposure. Mixing coefficients over a range of two 
orders of magnitude correspond to a large spread of neutron exposure. For larger 
mixing efficiencies the pocket containing s-process enriched material becomes 
broader. However, the mixing coefficients that reproduce neutron exposures 
estimated in are much smaller than what is predicted by the full stellar 
evolution model with rotation. We attempt to attribute this to unaccounted 
processes of angular momentum transport like magnetic fields (see Heger etal., 
these proceedings) in the AGB models. The simple parametric models indicate 
a possible solution for the apparent spread in s-process efficiency reported from 
pre-solar grains (Lugaro etal., 2002a) as well as from AGB and post- AGB stellar 
observations (Van Winckel & Reyniers, 2000; Busso etal., 2001). 
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